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By Stanley U. Bens cot er and Richard H. MacNeal 


SlMdARY 


Using the Cal-Tech analog con 5 >uter, structural analyses have heen 
made for four straight multicell wings. Wings with aspect ratios of 2 
and 4 with rectangular and biconvex cross sections have heen considered. 
The wings are siq)ported rigidly along two lines at the faces of the fuse- 
lage. Concentrated loads are applied at the intersection points of the 
rihs and spars. The effects of shearing strains in the ribs and spars 
are included. Deflections and all internal force quantities have been 
recorded as well as vibration modes and frequencies. 


INTRODircTION 


The structural analysis of a thin multicell wing of low aspect ratio 
presents a rather difficult problem if one wishes to obtain better accu- 
racy than can be obtained from elementary beam theory. Foiir such wings 
have been analyzed on the Cal-Tech analog computer for various static 
loads and in vibrational motion. Wings with rectangular and biconvex 
cross sections with aspect ratios of 2 and 4 have been analyzed. The 
wings are ass\imed to extend through the fuselage and to be supported 
rigidly along two lines at the faces of the fiaselage. 

VHien cross sections of a thin wing have a horizontal axis of sym- 
metry, the wing deforms \inder load in the manner of a plate. Analogous 
circuits for elastic plates were given in reference 1. The stmctural 
theory and analogo\is circuits for multicell wings are given in refer- 
ence 2. The present paper is devoted to a presentation of the results 
of con^jutations based i:ipon the method of reference 2. Con 5 >arisons with 
elementary beam theory are given wherever possible. In the case of 
bending loads the variation from beam theory becomes appreciable only at 
very low aspect ratios. In the case of torsional loads the variation 
from beam theory is more pronounced because of the occurrence of normal 
stresses due to warping restraint. The graphical illustrations represent 
only a small portion of the data which were obtained from the conputer. 
The con^Jlete resoilts of the computing work have been omitted for publica- 
tion but are given in the tables in the manuscript copy of this report. 
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This copy is available for loan or reference in the Division of Eesearch 
Information, National Advisory Committee for Aeronautics, Washington, D. C. 
The present investigation was conducted at the California Institute of 
Technology -under the sponsorship and -with the financial assistance of the 
National Advisory Committee for Aeronautics. 


SYMBOLS 


A 

A* 

Ain 

c 

d 

Pmn 

E 

f 

G 

I 

Iqi 

J 


area of a cell 

transformed area used in shear-flow calculations 
area of shear -web of ith rib in nth bay 
area of shear weh of jth spar in mth bay 
chord of wing 

transfonnation or carry-over factor 
bending stiffness of ith rih at jth spar 

bending stiffness of jth spar at ith rib 
twisting stiffness of a panel 
Young's modulus 
frequency 

shearing modtilus of elasticity 
depth of jth spar 

moment of inertia of skin per unit of -width 
manent of Inertia of jth spar 
total moment of Inertia of cross section 
torsion constant 


L 


distance from support line to -wing tip 
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M 

Mlj 

Mji 

P 

PiJ 

<1 

(^j)l 


"bending moment 

"bending moment in ith rl"b at jtb. spar 
"bending moment in Jth spar at ith rib 
concentrated load 
load at ith rlh and Jth spar 


shear flow 

shear flow in weh of jth spar 
statically detennlnate part of 




indeterminate part of 




^n 


R. 


R. 


R 


■xy 


^yx 

T 

Tb 

V 

Vin 

Vjm 


cell ul a r shear flow in nth cell 
flexibility of spar 


flexibility of rib 
interaction flexibility, 

interaction flexibility, 
twisting moment 

twisting moment coai5)uted from beam theory 


transformer numbers 
shear 

shear in ith rib at nth bay 
shear in jth spar 


shear in jth spar at mth bay 
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V 

^•b 

Vy 

w* 

a 

^in 
P jm 

\a 


% 


statically determinate part of Vj 
deflection 

deflection from "beam theory 
deflection at ith rib and jth spar 
deflection due to shearing strains only 

deflection of a beam from difference equations 
ratio of length, to width of a wall segment 
ratio of depth to width, of weh of jth spar 
rotation of normal in ith rih in nth hay 
rotation of normal in jth spar in mth hay 
jtmp in a function across ith rih 
Jun^) in a fimctxon across jth spar 
width of structure associated with ith rih 
width of struct\ire eissociated with jth spar 
width of mth hay between ribs 
width of nth hay between spars 
Poisson's ratio 

eq.\iivalent chordwise Poisson's ratio defined by 
eq.iaation (l) 

equivalent spanwise Poisson's ratio defined by 
equation (l) 

normal stress 

spanwise normal stress 

normal stress con5)uted from beam theory 
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2^ stmmiation about nth cell 

T shearing stress 

Tjj shearing stress cong)uted from beam theory 


DESCEIPTION OF SIRUCTURIS 


In .order to arrive at a reasonable set of dimensions for the struc- 
tures to be analyzed a few design con 5 )utations were made by using ele- 
mentary beam theory. The complete missile was assumed to weigh 20 kips 
and to be designed for an -ultimate load factor of 4. The allowable nor- 
mal stress in the skin -was assumed to be 20 ksi and the allowable shearing 
stress in the spar webs -was assumed to be 12 ksi. The loading -was assumed 
to be uniformly distributed. The thickness of the rib webs was arbitrarily 
ass\jmed to be one-half of the value determined for the spar webs. The 
resulting structixre for an aspect ratio of 2 with a rectangular cross sec- 
tion is shown in figure 1. 

In order that the skin might carry a normal stress of 20 ksi with- 
out buckling it would be necessary to reduce the spacing of the spars to 
approximately one -half of that shown in figure 1. This woxald provide 
13 spars rather than 7 spars . The limitation to seven spars -was dic- 
tated by the amount of electrical, equipment which available. Each 
spar which is shown in figure 1 should be considered as being equivalent 
to two spars in the staructure as it would be bviilt. Connection angles 
and shear web flanges have been omitted for convenience. 

Details of the -wing -with a rectangular cross section and an aspect 
ratio of are shown in figiire 2. The structural chord of 72 inches has 
been retained for all of the -wings. Also the fuselage bay has been 
assumed to be one-third of the total structural span for all wings . The 
allowable shearing stress in the spar webs was reduced to one -half of the 
previous value or 6 ksi. In order to develop this strength in the deeper 
web some form of stiffening against shear buckling would be req-uired. 

The dimensions of the wings -with biconvex sections have been chosen 
more or less arbitrarily to give about the same strength as that provided 
by the rectangialar cross sections. The dimensions of the -wing of aspect 
ratio 2 are shown in figure 3* For the wing having an aspect ratio of it 
and a biconvex section the dimensions are shown in figure it. 

In order to record the confuted results in -tabular form a numbering 
system for points on the plan form has been adopted as sho-wn in figure 5- 
A number has been assigned to each point at which some qioantity is measured. 
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In presenting the res\ilts of the analyses it is only necessary to give 
the results for one quadrant of the plan form since the layout of the 
structure has double symmetry in its plan. The points at which concen- 
trated loads are applied are shown in figure 6 . The deflections and 
tending moments are also measured at these points. The points at which 
shears are determined are shown in figure 7 * figure 8 the points 
are shown at which twisting moments are determined. These twisting 
moments may he considered as average values over the panels. 

A few sample calculations of structural constants are included to 
show in detail how these quantities are detennlned. The calculation of 
stiffness constants is made according to the formulas contained in ref- 
erence 2 . The wing is assumed to he constructed of an aluminm alloy 
having the material properties shown in table 1 . Calculations for the 
spanwise and chordwise hendlng stiffnesses at an interior point of the 
wing with rectangular cross section and aspect ratio 2 are given helow. 

For the spar; 



= (161.8 + 6,0) X 10^ 
= 167.8 X 10 ^ Ih-in.^ 


/l 6 l. 8 \ 

= I X 0.3 = 0.289 

V167.8/ 
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For the rlh: 



^ ^ (3-^r lO.if X 10^ 1 , ^ 

D. . = 16 X 0.16 X — X + — X 0.07 X (3.68)^ X 10.4 X 10° 

2 0.91 12 


= (216 + 3) X 10^ 


= 219 X 10^ Ih-in.^ 


= 



X 0.3 


0.296 


It is cortvenlent in the con5>utation of structural constants to intro- 
duce spanwlse and chordwlse Poisson ratios according to the following 
definition; 


Aj Ai 1 - [ 1 ^ 


The relation between bending moments and curvatures may be expressed in 
the following form: 


^4 A.p, 

T. i im j in 


^ ^ - -MlJ 


^i 


'J 


(2a) 


(2b) 
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In order to design the analogous circuits, equations (2) must he 
inverted to obtain formulas for the ciorvatures in terms of the bending 
moments. The coefficients of such equations may be referred to as 
flexibilities . Formulas for the flexibilities become : 


Rx = 


^i 


Dji(i - ^^ij^^Ji) 


















Dij(l - 


A san5)le calculation of the numerical values for the flexibilities at 
an interior point of a rectangular section with aspect ratio 2 are : 


1 - HijUji = 1 - (0.289) (0.296) = 0.914 




16 


167.8 X 10^ X 0.914 


= 0.1045 X (10)'^ 




12 


219 X 10^ X 0.914 


= 0.0600 X (10)"^ 
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12x0.296 , ,_6 

Rxy = = 0.0231 X (10) 

167.8 X 10^ X 0 . 9 li»- 


16 X 0.289 , , f , 

Ry^ = = 0.0231 X (10) ° 

219 X 10^ X o.9iij- 

A saii5)le calculation of the shearing constant of an interior spar is: 


16 X ^ X 10 ^ X O.llt X 3.68 


= 0.0303 X (10)"^ 


A calculation of the shearing constant of an interior rih is: 


^nGA-ia 12 x X 10 ^ X O.OT X 3.68 


= 0.0809 X ( 10 )"^ 


A sample calculation of the twisting constant for the skin is as follows; 


Gi = k X 10 ^ X 0.16 X = 4.72 X 10 ^ 


12 


\i _ 1 

IWi 16 X k.'J 2 X 10^ 


= 0.1589 X (10) 


-6 
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LOADING CONDITIONS 


Since the plan form of each wing has two axes of symmetry the 
loading may he divided into four types according to the symmetry con- 
ditions . Loads which are symmetrical about the midchord line are 
referred to as bending loads and loads which are antisymmetrlcal about 
the midchord line provide the torsional loading . Since the loads may 
be symmetric or antisymmetric about the plane of symmetry of the air- 
craft the following foijr types of loading may be considered; Symmetric 
bending, symmetric torsion, antisymmetric bending, and antisymmetric 
torsion. 

For the static loading conditions the loads are assumed to be 
applied as concentrated forces at the points of intersection of the ribs 
and spars. The loada are applied as a group of four forces, one in each 
quadrant of the plan fonn, to form a doubly symmetrical arrengement . In 
designing the analogous circuit and in recording the results of the analy- 
sis it is only necessary to consider one quadrant of the structure. A 
circuit was designed to correspond to the first quadrant of the plan form. 
Each point of this quadrant was loaded independently. The foin* types of 
loeiding were obtained by using the appropriate boundary conditions for 
the quadrant along the lines of symmetry of the plan form. 

Four wings, each to be loaded with loads having foiar types of sym- 
metry, provide 16 cases. In each case there are I6 loadings. Such a 
complete program of conputation woxiLd be too extensive to be jiostif iable . 
Hence it was decided to load all points in li- cases and to load only the 
tip points in the remaining 12 cases. The choice of loading points for 
the vario^as cases is shown in table 2 . In order to make a further saving 
in the calculation work for those cases wherein all points were loaded, 
the internal force quantities were not recorded when the foiir points in 
the fuselage bay were loaided. Deflections were, however, recorded in all 
cases . 


BOUNDARY CONDITIONS 


Along the leading edge of the wing the chordwise shear, chordwise 
bending moment, and chordwise twisting moment must vanish. Along the 
wing tip the spanwise shear, spanwise bending moment, and spanwise 
twisting moment must vanish. Along the line of support the vertical 
deflection must be zero. It has also been ass^lmed that the wing is 
damped along the support in such a manner as to prevent chordwise rota- 
tion of the normals to the elastic axis of the rib which lies over the 
support. This clamping, however, does not prevent the rib from warping 
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out of its own plane. This assuiiption that the normals in the rib do 
not rotate provides a relation between spanwlse chordwise bending 
moments over the si5)port. Since is zero at all points along the 

support the quantity may be set equal to zero in equations (2) . 

This gives the following formula for the chordwise bending moments in 
terms of the spanwlse bending moments: 


= ^^id 




(5) 


The center lines of the plan form are boundary lines for the quad- 
rant of the wing. Along these lines the boundary conditions are symmetry 
conditions and are determined from the known symmetry of .the applied load. 
It can be shown and it is obvious from physical intuition that the sym- 
metry of the deflection surface is the same as the symmetry of the applied 
load. From the known symmetry of the deflection surface the symmetry con- 
ditions of the various derivatives can be determined. From this informa- 
tion the symmetry conditions for the internal forces may be determined. 
This determination of the symmetry conditions for the internal forces 
follows the arguments, which are customary in elastic plate analysis. 


ANALOGOUS ELECTRICAL CIRCUITS 


Since all of the four wings to be analyzed have the same general 
arrangement of ribs and spars, the same arrangement of circuit elements 
may be \ised for all of the wings. The con5)lete circuit consists of three 
planar circuits which are shown in figvtres 9 to 11 . These circuits are 
similar to those given in reference 1 except that they Include the effect 
of shearing strains in the ribs and spars. The design of the circuits is 
in exact agreement with the methods of reference 2 . 

For each type of loading, changes along the planes of symmetry are 
required to satisfy the conditions of symmetry. These changes are indi- 
cated on the drawings of the circuits. Along the leading edge and wing 
tip the natural boundary conditions are automatically satisfied by the 
circuit. Along the s-upport the deflections and chordwise rotations of 
normals are made to vanish by groxaading the proper nodal points in fig- 
ures 9 and 11 . 
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AMLYSIS OF WINGS BY BEAM THEORY 


In presenting the results of the analysis it is of particular 
interest to conpare results obtained from the analog con 5 )uter with solu- 
tions obtained hy applying elementary beam theory to the wing. In tor- 
sion the ejcpression "elementary beam theory" is used to indicate that 
stresses and displacements at a cross section are computed from the sec- 
tion torque by using formulas from the solution of St. Venant. In 
bending action "elementary beam theory" is based on the assun^jtion that 
plane sections remain plane and in addition that cross sections are pre- 
served by ribs which are rigid in their own plane. This differs from 
the bending theory of St. Venant which predicts a chordwlse camber. All 
con 5 )arisons in this paper will be made with elementary beam theory. 

The shear flows in torsion were confuted by using the numerical pro- 
cedinre of reference 3* An example of this procedure is given in table 3 
for the biconvex section of the wing of aspect ratio 2. The solution is 
obtained by a process of successive corrections and checked by a cycle 
of iteration. Because of the symmetry of the cross section about a ver- 
tical axis it is only necessary to write the confutation for one-half of 
the section. At each step in the calculations it must be remembered that 
the shear flows in torsion are symmetrical about the center line. The 
torsion constant is also computed in table 3* torsion constants for 

all of the wings are given in table 1. 

QJhe shear flows in bending may be divided into a statically deter- 
minate part and an indeterminate part. This division may be indicated 
for the shear web of the jth spar as 


' (lj)s + (^})± 


w 


It is convenient to define the statically determinate shear flows in the 
top and bottom skin as being zero. The statically determinate shear flows 
in the spar webs aire confuted from shears which are proportional to the 
moments of inertia of the spars as follows; 



(5) 


The statically determinate shear flows are in equllibrlimi with the 
normal stresses anfl are also in equilibrium with the external loads. The 
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indeterminate shear flows are coii5)letely self-eq.uilihrating and are 
determined hy the condition of continuity of warping displacements. The 
Indeterminate shear flows consist of ce.llu3.ar shear flows which are anti- 
synmetric about the vertical axis of symmetry of the cross section. The 
indeterminate shear flow in a spar weh may he expressed as the difference 
between the cel lular shear flows in the two adjacent ce3Jn as follows: 


(ij)i “ * %-i. 


In writing equation (6) it is assumed that the jth spar lies between 
ceDJ. niimber n and cell number n - 1. 

The condition of continiiity of warping displacements leads to a 
difference equation governing the cell.ular shear flows as follows: 


-^'j+l'lnH + (^n°^)^n ' “-jV! = -^j+l(‘lj+l)s + ‘=^j(‘lj)s (T) 


Equation (7) is written for the nth cell. The coefficients are computed 
8LS ratios of length to width of the wall segments. The system of equa- 
tions corresponding to equation (7) has a matrix of coefficients which 
is the same as for the torsion case. Consequently the numerical procedure 
of reference 5 inay again be employed. An example of the calculation of 
shear flows in bending is given in table 4. Moments of inertia of the 
cross sections and of the individual spars are given in table 1. 


SOUECES OF ERROR 


Errors in the computed results due to inadequacy of the structural 
theory cannot be estimated. However, other soxnrces of error have been 
estinsted rovighly as follows: 

Error in conputer, percent ±1.0 

Error in meter, percent ±1.0 

Error in reading meter, percent ±0.5 

Error in conputing stiffness constants, percent ±1.0 

Error from oising difference equations (see figs. 10 and 3l1) 


An inspection of the data obtained from the conputer stiowed that the 
total bending moments and shears on various cross sections violated statics 
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"by amounts ranging from 1 percent to 3 percent. A factor to correct for 
this error was introduced at the same time that the electrical quantities 
were converted to structural quantities. The error in the conputer which 
is shown in the above table is an estimation of that error which remains 
after the correction for statics has been introduced. The error in the 
computer is due to parasitic effects in the con5>onents of the circuit. 

Some error is contained in the meter which is used for measiaring the 
volteiges or currents . Some additional human error is also Involved in 
reading the meter. The errors f;or the con5)uter and the meter, as given 
above, are applicable to the largest numbers which are given in the tables 
of data. For the smaller qiiantities the percentage errors may be much 
larger, even though the absolute size of the errors will be less. 

In the case of the rectangular cross sections there are no errors 
in the stiffness constants. However, in the case of the biconvex sec- 
tions there is a small error. A more significant error arises from the 
use of difference equations in expressing the structural theory. No 
specific value has been given for this error since there is a large dif- 
ference between the symmetric and antisymmetric cases. (See figs. 12 
anH 13.) In order to obtain an estimate of the error in deflections a 
ttniform beam was analyzed by solving difference equations and by solving 
differential equations . The central span and external spans were each 
divided into three parts in forming the difference equations . The beam 
was assumed to be loaded with a concentrated force at each tip. 

For the case of symmetric loads the ratio of deflections from dif- 
ference equations to deflections from differential equations is shown 
in figiare 12(b) . Here the error is seen to vary between 2 and 3 percent. 
The antisymmetric case is illustrated in figure 13(b) . In this case the 
error is more than five times as great as in the symmetric case. It is 
obvious that an accurate treatment of deflections in the antisymmetric 
case requires more stations along the span in forming the difference 
equations or a more accurate formulation of the structiiral theory in 
terms of difference equations . The effect of the use of difference 
equations on the distribution of internal forces in the multicell wings 
is not known. However, it can be shown that a3JL relations of statics 
between the internal forces and applied loads are satisfied. 

It is also of interest to estimate the error which would be made if 
the deflections due to shearing strains were neglected. This requires the 
use of specific cross sections. Such an estimate was computed for the 
wings with rectangular cross sections by using elementary beam theory. 

The ratio of deflections due only to shearing strains to the total deflec- 
tions, using differential equations, was conputed and plotted. For sym- 
metric loading the ratio is shown in figure 12(c) where it is seen that 
the error varies from 2 to 6 percent.. For the antisymmetric case, which 
is shown in figinre 13 (c), the ratio is approximately five times as large. 
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TREATMEOT OF DATA FROM COMPUTING MACHINE 


A sin5)le check which can he made on the results from the con?)uter 
consists of con?)aring the total shear, hendlrig moment, or torque on a 
cross section with the value given hy statics from the external loads. 

When this check was made errors were found which ranged from approxi- 
mately 1 to 5 percent. These errors were due to parasitic effects in the 
con^juter and not to the use of difference equations for the structural 
theory. The error did not show large variations over a given planar 
circxiit and hence it was decided that the voltages and currents in each 
planar circuit shoxild he conrected hy a single scalar factor. 

The correction factors were determined hy applying statics to the 
sections shown in figure l 4 . For hendlng loa^ the correction factor 
for hoth of the hending-moment circuits was determined to satisfy statics 
for the bending moment over the si^jport. The shear circuit was corrected 
to satisfy statics in the first hay outboard from the support. With tor- 
sional loads the hending-moment circuits were corrected to give the cor- 
rect section torque in the first hay outboard from the support. The 
shear circuit was corrected to give the correct total shear on section D-D 
as shown in figure l 4 (h) . Since the shears in the spars contribute to the 
section torque it is necessary to connect the shears first and the twisting 
moments suhsecLuently . 


DEFLECTIONS DUE TO BENDING LOADS 


The resiilts of the confutations on the analog computer are given in 
the tables available on loan from NACA. From these tables a small portion 
of the data has been taken to prepare illustrations of cases which seemed 
to be most Interesting. For all of the bending cases which are Illustrated 
the loading consists either of two comer loads at the tip or a single 
load on the central spar at the tip. 

The chordwise distribution of deflections at the tip is shown in fig- 
ure 15 for the wing of aspect ratio 2 with rectangular cross section. 

Both types of bending loads are applied symmetrically. For the central 
load the effect of shearing strains in the rib is clearly evident. A 
conf arisen is made with elementary beam theory using both differential 
equations and difference equations. The solution for a plate of infinite 
width with a uniform line load at the tip is also shown. This solution 

p 

is obtained from the solution for a beam by multiplying by 1 - p . It 
may be seen that the average deflection of the wing structure is somewhere 
between the values predicted by the two beam theories. 
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The chordwlse distributions of deflections for the wing of aspect 
ratio 2 anfl biconvex section are shown in figure 16. In this figure a 
centreil load is applied at the tip and a con5>arison is made between the 
symmetric and antisymmetric cases. Both cases show fairly good agree- 
ment with beam theory. 

Chordwise distributions of deflections along each of the ribs eire 
shown in figure I7 for both wings of biconvex section. The graphs are 
plotted as a ratio of con^nited deflection to the deflection from beam 
theory . 

The spanwise distributions of deflections for the leading-edge spar 
and the central spar are shown in figure I8 for the wing of aspect ratio 2 
with rectangular cross section. A coniparison is given between the sym- 
metric and antisymmetric cases. 


SHEAES IN SPARS DUE TO BENDING LOADS 


The chordwise distribirbions of shears in the spars at each outboard 
bay are shown in figure 19 for both wings of rectangiilar section. The 
ratio of shear from the conputer to shear obtained from beam theory is 
plotted. This ratio may also be regarded as a ratio of shearing stresses. 

For the wing of sispect ratio 2 with biconvex section the shears in 
the spars are illtistrated in figure 20. A coirparison is made of the dis- 
tribution for a central load and for comer loads. 

The spanwise distribution of shears in the leading-edge spar and 
the central spar are shown in figure 21 for a wing of aspect ratio 2 with 
rectangular cross section. Values obtained from beam theory are also 
shown. 


SHEARS IN RIBS DUE TO BENDING LOADS 


A single illustration of shears in the ribs is given in figure 22 
for the wing of aspect ratio 2 with rectangular section. The chordwise 
distribution of shear in all of the ribs is shown. These distributions 
are given for symmetrical cases of a central load and for comer loads. 


BMDING MOMENTS IN SPARS DUE TO BENDING LOADS 


From the spanwise bending moments one nay con^jute the spanwise nor- 
mal stresses. The ratio of spanwise normal stress from the conputer to 
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the same stress from "beam theory was confuted for all four wings for 
syrametric tending loads. The stresses were confuted for points over the 
support and are shown in figure 25 for a central tip load sind in fig- 
ure 24 for comer loads. It is of interest to note that the maximum 
variation from beam theory for the wing of aspect ratio 4 is only 3 per- 
cent while the variation for the wing of aspect ratio 2 is I5 percent. 

For each of the wings of aspect ratio 2 the normal stresses over 
the remaining rihs are shown in figvire 25 for a central load. In regions 
of smaller stress it is seen that the percentage variation from team 
theory becomes larger. 

The spanwlse variations of bending moment in the leading-edge spar 
and central speir are shown in figure 26 for the beam of aspect ratio 2 
with rectangular section. The effects of both comer loads and a central 
load are conpared with beam theory for the symmetric case. The corre- 
sponding antisymmetric case is shown in figure 27- 


BENDING MOMENTS IN BIBS DUE TO BENDING LOADS 


An illustration of chordwise bending moments in all of the ribs is 
shown in figure 28 for a wing of aspect ratio 2 with rectangular section. 
The effects of a central load and of comer loads are conpared for sym- 
metric cases. 


TWISTING MOMENTS DUE TO BENDING LOADS 


The chordwise distributions of twisting moments in the various bays 
are shown in figinre 29 for the wing of aspect ratio 2 with rectangular 
section. A conparison is made between the effects of a central load and 
comer loads. It may be seen that the twisting moments act in opposite 
directions for the two types of loads in all bays which are outboard from 
the fuseleige. 


DEFLECTIONS DUE TO TORSIONAL LOAD 


In preparing the illustrations for torsional loading the only loading 
which has been considered is that of a pair of equal and opposite forces 
at the comers . The chordwise distribution of deflections at- the tip is 
shown in flgtrre 30 for the wing of aspect ratio 2 with rectangular section. 
The symmetric and antisymmetric cases show very close agreement. The maxi- 
mum variation from beam theory is about 10 percent. 
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The chordwlse distrihutlon of deflections for symmetric loads is 
shown in figure 5I for the wing of aspect ratio 2 with a biconvex section. 
The curvature of the rih in this case is more noticeable than in the pre- 
vious figure. 

For the wing of aspect ratio it with biconvex section the chordwise 
distributions of deflections for synmetric loading are shown in flgture 52. • 
A conparlson with beam theory is given. The curvature of the tip rib is 
noticeable . 

Spanwise deflections of the leading-edge spar, and also the adjacent 
spar, are shown in figure 33 for file wing of aspect ratio 2 with a rec- 
tangulEir section. Symmetric and antisymmetric cases are conpared with 
beam theory. The same quantities are illustrated for the biconvex sec- 
tion in figure 3 ^- 


SHEARS IN SPARS DUE TO TORSIONAL LOAD 


The chordwise distributions of shears in the spars are shown in fig- 
ure 35 for both wings of aspect ratio 2 . Ratios of shearing stress as 
confuted to shearing stress from beam theory are shown. In almost all 
positions of the plan form the ratio is larger than unity. This is due 
to the torsion-bending action in which nornal stresses arise becavise of 
torsional loading. 

The spanwise distribution of shear in the spars is shown in fig- 
ure 36 for the wing of aspect ratio 2 with biconvex section. Large varia- 
tions from beam theory are evident. 


SHEARS IN RIBS DUE TO TORSIONAL LOADS 


The chordwise distributions of shears in the ribs are shown in fig- 
ure 37 for the two wings with rectangular section. A cougar Ison is shown 
for the two aspect-ratio cases for symmetric loads. 


BENDING MCMENTS IN SPARS DUE TO TORSIONAL LOAD 


For purposes of illustration the bending moments in the spars have 
been converted to normal stresses. The chordwise distributions of span- 
wise normal stress over the sioppoid; are shown in figure 38 for the fom: 
wings with symmetric loads. This lUiistration shows that the normal 
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stresses due to torsion are more than twice as large in the wings of 
lower aspect ratio. Similar results are given for the antisymmetric 
case in figure 39- 

The variation of spanwise normal stress along the leading-edge spar 
is shown in figure 40 for the wings of aspect ratio 2 . Results are shown 
for hoth types of cross sections with symmetric and antisymmetric loading. 
Corresponding stresses for the wings of aspect ratio 4 are shown in 
figure ifl. 


BENDING MOMENTS IN RIBS DUE TO TORSIONAL LOADS 


The Lending moments in the ribs have again been converted to normal 
stress. The distribution of chordwise normal stress in the various ribs 
is shown in figure h -2 for the wing of aspect ratio 2 with rectsingular 
section. S3raimetric loading is considered. 


TWISTING MOMENTS DUE TO TORSIONAL LOADS 


The chordwise distributions of twisting moments for various bays 
are shown in figure L3 for the wings of aspect ratio 2 . The quantity 
which is plotted is the ratio of computed twisting moment to the value 
of twisting moment given by beam theory. In most of the bays the ratio 
mity. This is due to the torsion-bending action of the 


is less t] 
structvir^ 


EFFECT OF SHEARING STRAINS IN RIBS AND SPARS 


Throughout all of the analyses on the analog con5)uter the effect of 
shearing strains was taken into account. These shearing strains affect 
not only the deflections but also the distribution of internal forces. 

A few calculations were made on the assun5>tion of an infinite shearing 
stiffness for the ribs and spans. The effect of shearing strains upon 
the distribution of shears in the spars is shown in figure for the 
wing of aspect ratio h with rectangular section. The change in shear 
due to the existence of finite shearing strains is generally small and 
is less than 10 percent in almost all cases. 

The effect of shearing strains upon the spanwise normal stress over 
the support is shown in figure for the wing of aspect ratio 4 with 
rectangular section. Symmetric bending and torsional loads are considered. 
The effect of shearing strains upon the normal stress distribution is some- 
what less than the effect \jpon shear distribution. 
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ANALYSIS OF A GKEDWORK OF RIBS AND SPARS 


In the analysis of a multlcell Beam the calculation work consists 
of two parts. In the first par*t the normal stresses are con^juted from 
the flexure formula and a statically determinate set of shear flows are 
comp\ited. In the second part of the analysis the indeterminate cellular 
shear flows are cong)uted. Consideration was given to the possihllity of 
developing a similar method of analysis for a multicell plate. In such 
a method it would he assumed that the skin is cut to prevent it from 
carrying any twisting moment. Shears and Bending moments would Be cal- 
culated for the gridwork of ribs and spars. A subsequent calculation 
would Be made to determine cellular shear flows required to restore con- 
tinuity of the skin. 

The above method of solution is Based on the assvm 5 )tion that the 
distribution of Bending moments in the gridwork will Be approximately 
the same as the distrlBirtion in the multicell shell. In order to check 
this assun 5 )tlon an analysis was made on the analog conputer of the wing 
of aspect ratio 2 as a gridwork. The results are con^jared with the 
analysis of the conjjlete wing in figure 46 for spanwlse normal stress 
over the support. From this figure it cein Be seen that with Bending 
loads the distribution from gridwork theory agrees fairly well with the 
correct distribution for the wing. However, under torsional loading the 
discrepancy is very large. Fiarther research study would Be required to 
develop this method of analysis. 


VIBRATION MODES 


In order to determine vibration modes and frequencies on the ansilog 
computer it is necessary to replace the resistors By inductors. It is 
assumed that the inertia forces are provided By an equivalent set of con- 
centrated masses . Each mass is represented in the analogous circuit By 
a condenser which is connected to ground. An illustration with detailed 
explanation of the vibration circuit for a Beam is given in reference 4-. 
In references 1, 2, and 4- the analogous circ\ilts have Been designed By 
connecting the condensers to the circtiit at the deflection points. In 
the present calculation work this technique has Been modified in order to 
obtain Better accioracy. 

It can Be shown in certain Beam analyses Based upon the use of dif- 
ference equations that Better accuracy can Be obtained if the equivalent 
concentrated masses are located at points midway Between the deflection 
points than if 'the nsnses ane located at the deflection points . This 
improvement in accuracy can also Be shown to Be true for a rectangular 
sinply supported plate. Consequently the wing modes have Been conputed 
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by considering concentrated masses to act at points on the ribs anfl spars 
midway between the intersection points of the ribs and spars. This divides 
the total mass into approximately twice as many small concentrated masses 
as woiold be obtained by assigning the masses to the deflection points. 

This greater subdivision of the total mass should also tend to improve 
the accuracy. 

Symmetric modes for the wing of aspect ratio 2 with rectangular cross' 
section eire illustrated in figure IfT by means of contour drawings. The 
first bending mode in figure shows very little camber. The mode 

shape is the same as would be obtained from beam theory. The second 
bending mode, as shown in figure 47(b), consists almost entirely of chord- 
wise bending. This camber mode cannot be predicted by beam theory. 

In the bending theory of Saint Venant the chordwise camber is obtained 
by multiplying the spanwlse curvature by Poisson's ratio. If this rela- 
tionship should tend to hold approximately true under various loading con- 
ditions it would provide a very convenient way of determining elastic cam- 
ber effects in aeroelastic analyses. A consideration of the first two 
bending modes as shown in figure 47 shows clearly that no relationship 
between spanwlse and chordwise curvatures can be assumed to hold vinder 
various loading conditions. 

The first torsion mode is shown in figure 47(c) . The frequency for 
this mode lies between the two values for the first two bending modes. 

The third bending mode is shown in figure 47(d) . 

The meaexired frequencies of the lowest modes in bending and torsion 
were con?)ared, for all four wings, with frequencies computed from beam 
theory without shearing strains. For the bending cases the variation 
from beam theory was less than 1 percent. The omission of shearing 
strains in the beam analysis compensated for the Poisson ratio effect in 
the wing. In the torsion cases the discrepancy was somewhat greater for 
the wing of aspect ratio 2. The wing of aspect ratio 4 showed a variation 
of 1 percent while the wing of aspect ratio 2 showed a frequency 8 percent 
higher than beam theory. 

For the wing of aspect ratio 4 with rectangular section two symmetric 
modes are shown in figinre 48. The first torsion mode is shown in fig- 
inre 48(a) . The first three bending imodes for this wing showed very little 
evidence of cairiber. A small amount of camber is seen in the third bending 
mode which is illustrated in figure 48(b) . 

The torsion modes can be lll\astrated more clearly by drawing deflec- 
tion diagrams for the various ribs and spars. The deflections of the 
spars in the first symmetric torsion miode are shown in figure 49 for the 
wing of aspect ratio 2 with rectangular section. The cinrves do not show 
the curvat'ure that would be obtained from beam theory. For the same mode 
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the deflections of the ribs are shown in figure 50. The curves show an 
appreciable amount of camber which is not predicted hy heam theory. 

The deflections of spars in the first symmetric torsion mode are 
shown in figure 5I for the wing of aspect ratio 4 with rectangular sec- 
tion. For the same mode the deflections of the rihs are illustrated 
in figure 52 . For the wing of aspect ratio If the deflections are more 
nearly like those which woiild be obtained from heam theory. 


CONCLUSIONS 


In preparing the illustrations for this paper it has only been pos- 
sible to consider a very small portion of the data which were conputed. - 
The reader will find that the data in the tables available on loan will 
permit him to make a more thorou^ study of any particular case in which 
he may be interested than could be given in the limited space herein 
available. However, a few general statements can be made. 

Under bending loads all of the wings show fairly good agreement with 
beam theory. The wings of aspect ratio 2 show poorer agreement than the 
wings of aspect ratio 4 but sti l l eigree, in general, more closely than one 
wo\ild be inclined to esjpect. The spanwise nojrmal stresses show about the 
same percentage variation from beam theory as the deflections . The shesiring 
stresses in the skin and spars, however, show much poorer agreement than 
the normal stresses. 

Under torsional lo6iding the w ing s show poorer agreement with beam 
theory than under bending loads. This is apparently due to the fact that 
torsional beam theory contains no effect of warping restraint (torsion- 
bending effect) . The warping restraint which actually exists causes nor- 
mal stresses in the spars and brings about a redistribution of the internal 
shear flows. It is clear that there is need for a practical theory of 
torsion bending for multicell beams. 

For the wing of aspect ratio 2 the second symmetric bending- vibration 
mri flp consists almost entirely of chordwise camber. For the wing of aspect 
ratio the first three symmetric bending modes contained only a small 
amount of chordwise camber. From the results of these analyses it is not 
unreasonable to believe that for wings of very low aspect ratio the effects 
of chordwise elastic camber should be included in aeroelastic analyses. 

In the torsional modes of vibration, due to the presence of warping 
restraint, the deflections of the spars showed appreciably less ctirvature 
near the tip than would be obtained from beam theory. 
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The data considered herein required the use of the analog computer 
for 4 weeks of normal working time. Consequently It may "be concluded that 
the analysis of such con^jlex structures as those herein considered is 
practicable with large automatic computers. 


California Institute of Technology, 

Pasadena, Calif., October 20, 1952. 
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OlABLE 1 

ELEMEINTAEY STRUCTURAL PROPERTBES 


(a) Material properties 



("b) Properties of total cross section 


Rectangular section 

Biconvex section 

Aspect ratio 

2 

k 

2 

k 

hjt , inA .... 

88.1^ 

386 

75.^ 

318 

J j jLu « • • • « 

308 

1,192 

268 



(c) Relative moments of inertia of spars 







table 2 


LOADED POINTS FOR VARIOUS CASES 



Rectangular 

section 

(a) 

Biconvex 

section 

(a) 

Aspect ratio 

2 

k 

2 

k 

Symmetric tending 

A 

T 

A 

T 

Symmetric torsion 

A 

T 

A 

T 

Antisymmetric bending 

T 

T 

T 

T 

Antisymmetric torsion 

T 

T 

T 

T 



all points loaxied; T, tip points loaded. 
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OlABLE 3 

SHEAR ELOTS IN TCQRSION FROM BEAM THEQRI 


jBiconvex section; aspect ratio, 2 ; 
J = 1514-. 0 X 2 = 268 in . 3 



1.62 

















Pi pi 
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TABLE 4 


SHEAR BLOWS DDE TO VERTICAL SHEAR EROM BEAM THEORY 


[Biconvex section; aspect ratio, 2; total 
section sliear, 2 





I8.81t 

16.20 

9.4o 

1.79 

• • • • 

0.512 

0.440 

0.255 

0.049 



4.3lt- 

4.02 

3.06 

1.46 

(Ij)s 

0.1180 

0.1093 

0.0854 

0.0556 

d 

a (qj)s . . . 

(0.i4t8) 

O.OlTit 

(0.1369) (0.1431) 
-0.0150 0.0156 

(0.1092) (0.1207) 
-0.0091 0.0101 

(0.0575) 

-0.0019 


(0.1478) 

0.0024 

-0.0004 

-0.0001 

(0.1431) (0.1369) 
0.0065 
0.0009 0.0005 

0.0002 0.0001 

(0.1207) (0.1092) 
0.0082 
0.0009 0.0008 

0.0001 0.0001 



0.0030 

-0.0004 

0.0079 
0.0011 0.0004 

0.0091 

0.0010 0.0010 




0.0031 

0.0079 

0.0092 


(q.j)i .... 

-0.0062 

-0.0048 

-0.0015 

0.0092 

•Ij 

0.1118 

0.1045 

0.0821 

0.0428 



0.485 

0.420 

0.251 

0.062 
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L)ne of 
suppoi-f 


("b) View A-A. 


OM’ 






0J4' 


36 ' 


45' 


(c) Section B-B. 

Figure 5 .- Wing witli Biconvex cross section. Aspect ratio, 2 
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Line, of support 
C Face of fuse/o^) 


Leocffn(^ ed^&of v\/}n^ 



2/^ (3/ 


Midchotd of win^ 
Plane of symmetry of oitcrafi 


Figure 5-- numbering of points at vbicb quantities are measured. 
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Figure 6.- Points at whlcli concentrated loads are applied nrifl deflections 

and 'bending moments are measured. 
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Figure 7.- Points at -which spanwise and chord-wise shears are determined. 


0 

6 

26 

46 

66 

86 

0 

4 

24 

44 

64 

B4 

c 


22 

4B 

62 

82 


Figure 8 .- Points at -which -twisting moments are determined. 


















Ramovti for 
iymm^frtc cose 


Shorf fhese nodes 
for torsion 


Figure 9-- Circuit for shears and deflections. 



Dhorr i-hese nodes 
for symmefric ca»c 


f^i f f f 9 

for forsion 


Figure 10.- Circuit for 8pamd.se ‘bending moments and chordwiee twisting 

moments . 




RejT^VG Rot 
symmetric, Cos& 


Remove, for 
forston 


Figure 11.- Circuit for chordwlse 'bending moments and spanwlse twisting 

moments . 
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P 




A/?.= 




(c) Deflections due to shear for rectangular cross section. 
Figure 12.- Deflections from beam theory for symmetrical loads. 
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(a) Loaded beam. 




A.R. 


A.R. 


(c) Deflections due to shear for rectangular cross section. 
Figure 15.- Deflections from beam theory for anti symmetrical loads. 
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load 

(a) Be n d in g load. Bending moment on section A-A Is 2PL. Shear on 

section B-B Is 2P. 



ioQd 

(h) Torsional load. Torque on section C-C is Pc. 

Shear on section D-D Is- P. 

Figure l4.- Control sections for correcting electrical quantities to 

satisfy statics. 
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(a) Loading points. P = 1 kip. 


• (b) Half cross section. 




(c) Deflections. 

Figure 15 •- Chordwlse deflections at tip. Aspect ratio, 2 . 

Symmetric bending. 




j 3 S 7 


Diffet^nctt mfuerfion^ 

m^u^Hons 


(c) Deflections. 

Figure l6.- (Jhordwise deflections at tip. Aspect ratio, 2. Central 

load at tip. 







(b) Relative values of deflections. 

Figure I7.- Chordwise deflections at ribs. Biconvex section. 

Symmetric bending. 









shears in spars. Symmetric "bending . 
section. 







(c) Antisymmetric loads. 


Figure 21.- Spanwise distribution of sheeirs in spars. Aspect ratio, 2. 

Rectangular section. 
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(c) Aspect ratio, i»-. 

Figure 23 -- Spanwise nonral stress over support. Symmetric bending. 

Central load.. 
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(a) Ijoading points. 


£x 


£x 



/.oz 

t.oo- 

.98- 




Redonquhr 

^ecNon^ j. 








Biconvex 

section 



(c) Aspect ratio, 4. 

Figure 2h.- Spanwlse normal stress over support. Symmetric "bending. 

Comer loads . 
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(a) Loading points. 




mb 3 







(Id ) Twisting moments . 

Figure 29.- Twisting moments on chordwise cross sections. Aspect ratio^ 2 
Symmetric "bending. Rectangular section. 
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(a) Loading points. P = 1 poimd. 


3^m theo 


Beam fhmory 


leo/r? theory 


(Id) Shear dlstrlhutlons . 

Spanwlse distribution of shear In spars. Aspect ratio, 2 . 
Symmetric torsion. Biconvex section. 


















1 pound. 
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(b ) Relative values of twisting moments . 

Figure k- 3 -- JTwisting moments on chordwlse cross-sections. Symmetric 

torsion. Aspect ratio, 2 . 
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(a) Loading points. 



(l) Relative values of sliearing stress. 


Figure 44.- Effect of sliearing strains upon distribution of shears in 
spars. Aspect ratio, ll-. Symmetric bending. Rectangular section. 
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(a) Loading points. P = 1 pound. 



g./ 1 

. ^ h/ifh sh&ar ^^frairr^ ^ ^ ^ 

'*9 / 

iV/fhout shear strain 

/ 3 5 7 

(b) Ssnnmetric bending. 









' 


t — — 

h/ifh sh&ar ^ 

'train 






h//fhout sheoi 

r strain 



(c) Symmetric torsion. 

Figure 45-- Effect of sbeaxing strains upon spemwise normal stress over 
support. Rectangulax section. Aspect ratio, 4. 










(a) First symmetric ■bending 
mode; f = 14-8.0 cps. 


(b) Second synmetrlc bending 
mode ; f = 26 l cps . 


Figure ' 14-7 . - "Vibration modes. Rectangular cross section. Aspect ratio, 2. 
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(b) Third symmetric bending mode; f = I85 cps. 

Figure lv 8 .- Vibration modes. Rectangular cross section. Aspect ratio, 4 . 




















